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We perform classical Monte Carlo simulations of ultracold neutron transport through an 
absorbing-reflecting mirror system in the Earth's gravitational field. We show that the underly- 
ing mixed phase space of regular skipping motion and random motion due to disorder scattering 
can be exploited to realize a velocity filter for ultracold neutrons. The range of velocities selected is 
controlled by geometric parameters of the wave guide. Possible applications include investigations 
of transport and scattering dynamics in confined systems. 
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I. INTRODUCTION 

Neutrons with an energy of around 100 neV or below 
are conventionally referred to as ultracold. One of their 
important properties is total reflection from many mate- 
rials allowing for storage of ultracold neutrons (UCN) for 
long times up to their lifetime in material traps or sealed 
vessels. UCN are conveniently produced in a research re- 
actor by extracting the low energy part of the Maxwellian 
spectrum of thermal neutrons and by directing them with 
neutron guides of different forms and sections [1, 2]. Well 
known applications of UCN include measurements of the 
neutron electric dipole moment [3, 4], the neutron half- 
life time [1, 2] as well as observation of the UCN quantum 
states in the Earth's gravitational field [5-8]. 

UCN can also serve as a powerful tool for study- 
ing quantum chaos [9], scattering of ultracold fermions 
[10, 11] and transport in confined geometries in the pres- 
ence of only short-ranged interactions while eliminating 
distortions by Coulomb correlations and interactions. A 
major complication is, however, the loss of neutrons due 
to absorption by the confining material. The loss oc- 
curs primarily if the neutron velocity perpendicular to 
the boundary surfaces (v±) exceeds a characteristic criti- 
cal velocity (v cr ) for total reflection of the confining mat- 
ter. Absorption is enhanced in wave guides with surface 
or bulk disorder as scattering at rough surfaces or im- 
purities can increase v±. The loss can be minimized by 
pre-selecting only such neutrons whose total velocity does 
not exceed the critical velocity of the waveguide walls. 
However, in the experiment without an additional selec- 
tor such a condition is hard to achieve since the velocity 
spectra in the initial UCN beam usually shows nonzero 
probability of neutrons with velocities much larger than 
v cr . In this paper we propose and simulate the efficiency 
of a band pass velocity filter for UCN. The proposed fil- 
ter consists of a neutron absorbing-reflecting mirror setup 
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similar to the one used for the investigation of neutron 
quantum states [5, 7]. We focus on system parameters for 
which the system approaches the classical limit and quan- 
tum effects can be neglected. For sampling the classi- 
cal phase space we employ the classical trajectory Monte 
Carlo (CTMC) method and derive analytic estimates for 
the velocity cut-offs of the filter. 



II. SET-UP OF THE FILTER 

Figure 1 shows a schematic picture of the proposed 
absorbing-reflecting mirror system (ARMS). The UCN 
beam from the neutron guide passing a collimator enters 
the ARMS which consist of two horizontal borosilicate 
glass plates, or mirrors, separated by a distance W. W 
is considered to be large enough (W > 160 fim) to sup- 
press quantum effects (such as diffraction and localiza- 
tion) and to allow for a classical transport description. 
The two glass plates are separated by two brass spacers 
which are rectangular blocks with height W and length 




FIG. 1. Schematic view of ARMS proposed as the veloc- 
ity filter consisting of a glass mirror (bottom) and a rough 
glass plate with surface disorder inducing large angle disorder 
scattering and absorption (top). The incident neutron beam 
leaves the neutron guide through a collimator. For the spacers 
position see Fig. 8(b). 
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FIG. 2. A typical surface of the rough upper mirror coated 
with Ni measured with the scanning electron microscope. 
(The measurement was performed at the USTEM facility 
[14].) The disorder has amplitude A ~ 3 \ivn and a corre- 
lation length l cor ~ 2 \±vn. 



equal to the mirror length l x (see Fig. 8(b), not shown 
in Fig. 1) placed symmetrically at the sides of the glass 
plates. The inner surface of the upper mirror is rough 
where the roughness is characterized by the amplitude A 
of the disorder and its correlation length l cor . A typical 
surface profile measured with a scanning electron micro- 
scope is shown in Fig. 2. The inner surface of the lower 
mirror is flat allowing for specular reflections. Only the 
neutrons which pass through this device without being 
absorbed by the glass plates (or by the spacers) are de- 
tected. The detector covers an area of W x l y (Fig. 1). 

To understand how this device can operate as a ve- 
locity filter, we analyze the classical phase space of this 
system with mixed regular- chaotic dynamics. A neutron 
traveling between two mirrors has 3 degrees of freedom 
and its Hamiltonian is given by 



H(p x ,Py,p z ,x,y,z) 



2m 2m 



PL 

2m 



+ mgz + V(x, y, z), 

(i) 

where g is the Earth's gravitational constant and 
V(x, y, z) describes the potential of the walls approxi- 
mated by: 



V(x,y,z) 



0, if < z < £(x,y) 
Vp, otherwise, 



(2) 



where £(x,y) describes the disorder landscape of the up- 
per disordered surface and Vp will be discussed later. 
The Hamiltonian in Eq. (1) governs the trajectories in 
6-dimensional phase space. A useful method for the 
visualization of trajectories lying on the 5-dimensional 
isoenergy surface in multi-dimensional phase space is the 
projected Poincare surface of section. We start from a 
four-dimensional surface of section given by a cut at the 
system exit x = l x with periodic boundary condition in x 




FIG. 3. The projection of the Poincare surface of section on 
the (p z , z) plane of the phase space for neutron transport con- 
sists of a regular island (red area at the center) representing 
skipping trajectories, and of an irregular part of scattering 
trajectories induced by the surface disorder with roughness 
amplitude A. The blue area with z ranging from W to W -\- A 
indicates the region where neutrons scatter off the upper sur- 
face. 



direction. We then project the surface of section onto the 
(Pzi z ) plane for p x > 0. Fig. 3 clearly displays two types 
of motion (note that due to time-reversal symmetry the 
surface of section for p x < is equivalent). The large 
regular island at the center corresponds to the neutrons 
repeatedly bouncing off the lower mirror in the Earth's 
gravitational field (skipping motion, [12]) and traveling 
through the filter without touching the upper rough wall. 
For this motion the degrees of freedoms are separable and 
the trajectory intersection lies on a smooth curve in the 
(Pz-, z) plane. The equation for this curve p z = p z (z) fol- 
lows from the Hamilton equation with the Hamiltonian 
H(p z ,z) = H z o being the integral of motion and inde- 
pendent of the other two degrees of freedom: 



Pz 



±m\J 2g(h 



(3) 



where h is a maximal height that a neutron can achieve 
for a given energy H z o — mgh. The region of skipping 
motion in phase space is in contrast to the region of 
chaotic motion of those neutrons that scatter off the dis- 
ordered upper surface creating the chaotic sea in phase 
space. In the latter case, due to the roughness, the de- 
grees of freedom are not separable and additional inte- 
grals of motion, besides the total energy, no longer exist. 
Therefore, the intersection of randomly scattered trajec- 
tory with the projected surface of section does not form 
smooth curve and shows chaotic behavior instead [15]. 

The interaction of neutrons with the surface of the 
glass walls as well as the spacers is described by a Fermi 
pseudopotential Vp which characterizes the low-energy 
nuclear scattering of surface atoms [1,2] and can be writ- 
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ten as 



III. MONTE CARLO SIMULATION 



\'i \'it <V, = 2 ^-N(b-ib i ), 



m 



where m is the neutron mass, N is the scattering center 
density, b is the coherent nuclear scattering length, and 
the imaginary part of the potential is determined by 



(4) To explore the implications of the phase space struc- 

ture analyzed above for the neutron transport through 
the ARMS, we perform a full three-dimensional Monte 
Carlo simulation using the geometric parameters of typ- 
ical set-ups and realistic initial distributions of the neu- 
tron beam. 



(5) 



where v is the neutron velocity and o\ is the absorp- 
tion scattering cross section due to the interaction of 
neutrons with nuclei of surface atoms. The absorption 
results from neutron capture by nuclei accompanied by 
emission of 7-rays or charged particles. This effect is 
accounted for within our description of elastic neutron 
transport by introducing the effective imaginary poten- 
tial Vj in (4). Within the classical simulation, it gives 
rise to a loss probability for the trajectories. 

The real part of the Fermi potential Vr determines the 
critical energy and velocity v cr of neutrons which are re- 
flected. If a neutron has a kinetic energy in the direction 
normal to the surface larger than the real part of the 
Fermi potential of the surface material (E± > Vr), it 
penetrates the wall and is considered to be lost. The sur- 
face roughness gives rise to nonspecular reflection and, 
thus, to transfer of longitudinal (Eu) to transverse {E±_ 
energy of neutrons. As a result, the total absorption 
will be enhanced by surface disorder. Consequently, neu- 
trons scattered off the upper wall, i.e. neutrons under- 
going chaotic scattering are effectively removed from the 
transmission. 

Even for E± < Vr the reflection probability from the 
mirror surfaces is below one, due to an exponential pen- 
etration of the wave function into the barrier ("tunnel- 
ing"). The reflection probability of the neutron from the 
step-like potential barrier follows as (e.g. see [2, 16]) 



R = 



k + k 



EY - (E ± - Vf) 1 / 2 

(E ± -v F y/\ 



pV2 



(6) 



where fco and k are the normal components of the wave 
vector in vacuum and in the medium. For the case Vj <C 
Vr (usually Vj is 10 4 times smaller than Vr) the reflection 
coefficient, to first order in Vj/Vr, becomes [2]: 



R= 1 - fi(E ± ) = 1 



E± 



v R y v r -e ± 



(7) 



with E±_ < Vr. The wall- loss probability per bounce is 
denoted by fi(E_\_). It should be noted that Vi /Vr = 
aiv/Anhb is independent of the neutron velocity since 07 
is proportional to 1/v [1]. Multiple scattering leads to 
preferential depletion of the chaotic region of the phase 
space since the number of bounces at walls for chaotic 
scattering trajectories is much larger than for skipping 
trajectories. 



A. Geometry of the set-up 



The two parallel absorber-mirrors are chosen to be 
square plates with l x = l y = 10 cm (see Fig. 1). To 
stay within the classical regime the distance between the 
two mirrors is selected to be W = 250 fim, in agree- 
ment with experiments (see, e.g. [5]). The mirrors 
are made of borosilicate glass and have a critical ve- 
locity v cr = 4.3 m/s and a loss factor (see Eq. (7)) of 
VijV R = 0.0015 [1]. The brass spacers (see Fig. 8(b)) 
have a critical velocity v cr = 5.3 m/s [1]. To simulate 
the landscape of surface disorder y) in Hamiltonian 
in Eq. (1) with realistic spatial correlations we use a con- 
volution method with a Gaussian correlation function: 



) (f (a, y)£(x + x', y + y')) = A 2 exp(- 



2/ 2 



)exp(- 



2/ 2 ' 



(8) 

The description of this method can be found in [17, 18]. 
The amplitude and the correlation length of the surface 
disorder in the simulations are chosen to be A = 3 \im 
and Iqov ~~ 

2 \im taken from the measured surface profile 
(Fig. 2). The distance between collimator and mirrors 
is Id = 18 cm, the distance between the lower collimator 
edge and the lower mirror surface is lh = 5.4 mm. The 
size of the collimator opening is y x z = 8cm x 2.7mm. We 
emphasize that the present simulation of the two-mirror 
setup taking into account also the transverse degree of 
freedom y is fully 3D unlike earlier two-dimensional mod- 
els for similar systems [7, 13]. 



B. Initial conditions 

The simulation requires the phase space distribution 
p(r, v, t = 0) of the UCN beam as it exits the collimator 
as initial conditions for the propagation. The most re- 
liable information about the distribution of the forward 
velocity component v x of the neutron beam at the exit 
of the neutron guide can be obtained from time-of-flight 
(TOF) measurements. In such a measurement, the neu- 
tron beam is periodically chopped (see, e.g. [19]). Af- 
ter this chopper, the neutrons traverse a neutron guide 
of known length before they are detected with the help 
of a neutron counter. The detector electronics also re- 
ceives one logical signal per turn of the chopper blade 
and thereby measures the time of flight of the neutrons. 
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FIG. 4. Velocity distribution f(v x ) measured by time-of- flight 
at the ILL facility for a straight, 2 m long cylindrical neutron 
guide with 4 cm radius. 
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The time-of-flight spectrum of the ultracold neutron 
beam at the Institute Laue-Langevin was measured with 
this method in 2008 [20]. The time-of-flight spectrum 
was corrected for the chopper offset. The opening func- 
tion of the chopper was deconvoluted and with the help 
of the distance-of-flight, the corresponding velocity spec- 
trum was extracted. The resulting normalized distribu- 
tion of the forward velocity in ^-direction (Fig. 4) will be 
used as initial distribution function. 

Although it is difficult to measure the transverse veloc- 
ity components v y and v z , part of the angular distribution 
P(0) of the beam are available (see, e.g. [1, 21]). P{0) 
depends on the ratio R of length to radius of the neu- 
tron guide. Since the angular distribution for the pipe 
which was used for the TOF measurement is unknown 
we perform our calculations for two distributions which 
correspond to the following limiting cases: in case of the 
ratio R ^> 1 the angular distribution is well approxi- 
mated by P(0) ~ cos(#) and in case of a short pipe with 
R > 1 the angular distribution is uniform. Because of 
the cylindrical symmetry of the pipe and the small an- 
gles involved, we consider the distributions P(0 y ) and 
P(0 Z ) to be identical and either cosine or uniformly dis- 
tributed. Since the inner surface of the cylindrical neu- 
tron guide is coated with Ni 58 , which has a critical ve- 
locity of 8.07 m/s, the velocity component normal to the 

guide surface at the pipe exit v± = * lv% + v\ should be 



less than 8.07ra/s. The initial vectorial velocity distribu- 
tion is randomly sampled from these distributions. The 
coordinate distribution of the neutron beam in the plane 
of the collimator opening is assumed to be uniform. 



C. Simulation results 

The classical equations of motion are integrated for a 
large number of initial conditions p(t = 0). Calculations 
start at the collimator opening with initial conditions 



FIG. 5. A sample distribution of 3000 trajectories undergoing 
chaotic scattering in the top view of the absorbing-reflecting 
mirror system (ARMS). The trajectories enter with initial 
conditions y — and v y = 0. Color scale in arbitrary units. 



randomly selected for each particle from the distributions 
discussed in section IIIB. The number of neutrons that 
reach the mirror slit entrance is kept fixed at 5 • 10 5 for 
the results presented in this section. The propagation in- 
side the slit includes absorption by walls and by spacers 
when the kinetic energy perpendicular to the surface is 
larger than the Fermi potential. The wall-loss probability 
per bounce Eq. (7) is included as a stochastic process for 
each collision with the wall. An example of the trajectory 
distribution between two mirror plates in top view of the 
system is shown in Fig. 5. The trajectories depicted here 
are calculated for 3000 different initial conditions lying 
inside the region of chaotic scattering. For better visual- 
ization we restricted here the initial conditions to y = 
and v y = 0. We also do not show the skipping trajecto- 
ries as they would coalesce to a single line (y = 0, in this 
case) and would decrease the contrast of the irregular 
scattering pattern. 

Results of the classical Monte Carlo simulations of the 
full 3D filter system are presented in Fig. 6, where distri- 
butions f(vi) of the absolute value of the total velocity 
(a) as well as the v x , v y and v z velocity components (c- 
d) at the detector position are shown for an initial cosine 
angular distribution. If a uniform initial angular distri- 
bution in the beam is used instead, the results hardly 
differ from the ones presented in Fig. 6. This is a con- 
sequence of the fact that the initial maximal angles are 
rather small (6 y ~ 26° and Z ~ 7°). Therefore, the uni- 
form and cosine distributions do not significantly differ 
from each other. 

The Poincare surface of section x = l x in the phase 
space projected on the (v z ,z) plane for the case of a 
cosine angular distributions (Fig. 7(a)) prominently dis- 
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plays a regular island of skipping motion with dense non- 
uniform distribution. The distribution of the chaotic 
scattering trajectories is barely visible in this plot since 
they spread out on a much larger region in phase space. 



D. Analytic estimate of the velocity band pass 

Each velocity component distribution can be decom- 
posed into contributions from skipping and chaotic scat- 
tering trajectories (blue dotted and green dashed lines in 
Fig. 6). For the skipping part of phase space it is possible 
to provide an analytical estimate for the velocity range, 
i.e. the minimum and maximum velocity of neutrons that 
can pass the filter. The shape of the distributions is more 
difficult to estimate since it strongly depends on the ini- 
tial velocities in the neutron beam. However, the velocity 
limits are independent of the initial distribution (pro- 
vided it lies within the skipping region of phase space). 

The band pass in forward direction, v x , can be cal- 
culated as follows. The maximum (minimum) velocity 
corresponds to the minimum (maximum) time of flight 
between the collimator and the glass mirrors 



— Id/ • 



(9) 



Moreover, neutrons have to enter the mirror system 
in the skipping region: neutrons do not touch the upper 
rough surface if the vertical velocity component at the 
collimator lies in the interval 



^2g8~z < v z0 < V / 2#TWj, 



(10) 



where 5z is the distance between the lower glass plate and 
an initial z-coordinate of the neutron at the collimator 
opening. Taking into account the motion in z direction 
allows to calculate the neutron time of flight between the 
collimator and the glass mirrors (Eg. (9), 



t± = (v z0 - v za )/g, 



(11) 



where v z q — v z (t — 0) is the initial vertical velocity and 
v zs is the velocity component at the mirror entrance slit, 



±\ v* -2g(5z + z s i). 



(12) 



Here z s \ is the vertical distance from the point where the 
trajectory enters the ARMS to the lower mirror. Eq. (9) 
and Eq. (11) gives v x as a function of v z q and Sz 



gU 



z0 ± \A 



v 2 z0 - 2g(5z + z s i) 



(13) 



For the fastest (slowest) neutrons the denominator in 
Eq. (13) has the smallest (largest) value corresponding to 
the -(+) sign. Also, the absolute value of v zs in Eq. (12) 
should be maximal which is reached for z s \ = and v z q at 
its maximum value obtained from Eq. (10). In addition 



Sz takes on the smallest (largest) possible value. The two 
trajectories with maximal and minimal v x are illustrated 
in Fig. 8(a) (side view) and correspond to the velocities 



gld 



W)-y/2gW 



\/2g(5z d 

o 



W) + y/2gW 



(14) 



where Sz up (Szdown) is the distance between the upper 
(lower) edge of the collimator opening and the lower glass 
plate (see Fig. 8(a)). Eq. (14) represents the analytic 
estimate for the velocity band pass for v x . Obviously, the 
velocity range in Eq. (14) is controlled only by geometric 
parameters of the system and can be tuned accordingly. 

The range of transmitted vertical velocity component 
v z for trajectories undergoing skipping motion can be ob- 
tained from the requirement that neutrons do not touch 
the upper wall inside the mirror slit, i.e., 



(15) 



The range of the v y component can be found from 
Fig. 8(b) representing the top view of the system. Here 
the trajectories with the largest angles that can still be 
detected are shown. In the presence of spacers we find 



\Vy\ < V x 



ly ~^ ly 
: 2L 



(16a) 



v yc is the size of the collimator opening in y- 



where l y 

direction. Without spacers the allowed velocity range 
is reduced to 



\Vy\ < Vx 



ly ~\~ lyc 

''2(i d + i x y 



(16b) 



since neutrons that would have been reflected from the 
spacers and stay inside the system are now lost. The v y 
range depends also only on the geometric parameters of 
the system and can be tuned as well. 

The limits Eq. (14 - 16) are depicted in Fig. 6(b-d) 
as dashed vertical lines and coincide with limits for the 
distributions of the skipping states (blue dotted lines) 
obtained from Monte Carlo simulations. It should be 
noted that the band pass filter limits for the v x and v z 
velocity components are separable, i.e. independent of the 
other Cartesian components of the velocity vector and 
only controlled by geometry. Only the v y limit depends, 
in addition to geometric parameters, on the maximum 
velocity in x direction. Consequently, the ARMS can 
serve as a band pass filter for the full vectorial velocity 
v, however only for skipping motion. 

It is now crucial to inquire the relative importance of 
skipping compared to chaotic trajectories and gather pa- 
rameters to enhance the relative weight of the skipping 
region of phase space for which the filter is operational. 
For the case depicted in Fig. 6, the fraction of chaotic 
scattering trajectories 



(17) 



6 




FIG. 6. Monte Carlo simulations for the components of the vectorial velocity distribution after passing the ARMS with 
parameters discussed in section III A. (a) Distribution of the absolute value of the total velocity f(v)\ (b-d) distributions of 
the components f(v x ), f(v y ) and f(v z ). The total distributions (solid red lines) can be decomposed into gravitational skipping 
(dotted blue line) and chaotic scattering (dashed green line) trajectories. Vertical lines correspond to analytical results for the 
band pass of the velocity components Eq. (14-16). The parameters of the filter: W = 250 fim, l x = l y = 10 era, A = 3 /xra, 
Iqot — ^ jJLinrLm 



where T CjS denotes the volume in phase space pertaining 
to either chaotic (c) or skipping (s) motion, is F c = 0.27. 
This fraction can, however, be tuned (see Fig. 9) by ex- 
ploiting the fact that chaotic trajectories executing many 
random bounces at the mirror are subject to enhanced 
absorption (see Eq. (7)). For example, for the distri- 
butions of Fig. 6, the absorption by wall collisions sup- 
presses the flux of chaotic scattering trajectories at the 
detector by 12%. In contrast, the flux of neutrons un- 
dergoing skipping motion is reduced only by 0.2%, be- 
cause of the drastically different number of collisions for 
chaotic scattering and skipping motion (for this exam- 
ple ~ 500 and ~ 10, respectively). By further increas- 



ing the number of bounces the phase space region as- 
sociated with chaotic scattering can be even more sup- 
pressed. Trajectories belonging to T c give rise to a broad 
velocity distribution at the detector extending well be- 
yond v cr and the band pass (Eq. (14-16), see Fig. 6). 
The low- velocity portion corresponds to trajectories with 
a large number of bounces. Since the total velocity of 
these neutrons is less than v cr they do not penetrate the 
walls. The leading loss mechanism is given by Eq. (7) 
which reduces the number of these neutrons but does 
not completely remove them. A simple estimate for re- 
flection losses (without taking into account the angular 
dependencies) is ~ exp(— N Vl /y R ), where N is number of 
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v z (m/s) v z (m/s) 



FIG. 7. Projection of the Poincare section x = l x on the (v z , z) plane for the filter system with (a) l x = 0.1 m and (b) 
l x = 0.2 m. Only the regular island is visible since the trajectories in the chaotic sea are spread out over a large volume in the 
phase space resulting in a much reduced density in the projected surface of section compared to the regular island. 



bounces and Vi/Vr = 0.0015. Strongly suppressing the 
amount of these low-velocity neutrons would require ~ 
2000 and more wall collisions. The high velocity part of 
the spectrum corresponds to neutrons undergoing only a 
few collisions with the rough wall. The normal velocity, 
v±, remains still less than v cr allowing these particles to 



side view 




FIG. 8. Trajectories representing the limiting velocities for 
transmission, (a) the v x component, the solid red line corre- 
sponds to the trajectory with minimal forward velocity, the 
dashed blue line to the one with maximal v x . (b) the v y 
component, the solid red line corresponds to the trajectory 
launching with the largest angle y still reflected by the spac- 
ers (depicted as black horizontal rectangles), the dashed line 
displays the same in case without spacers. 



travel through the mirror system without being lost. 

One more aspect to be discussed is the influence of 
spacers. Our Monte Carlo simulations indicate that 34% 
of the total number of trajectories (24% from skipping 
motion and 10% from scattering motion) are reflected 
from the spacers. We consider spacers to be flat and to 
be perfect specular reflectors, i.e. only v y changes sign af- 
ter scattering. Realistically the spacers do have a rough 
surface as well as other surface imperfections and thus 
may also give rise to randomization of the reflection an- 
gle. This can lead to an additional contribution to the 
irregular part of phase space. For example, skipping tra- 
jectories can become chaotic scattering trajectories after 
a single reflection from a spacer. Thus, in the foolowing 
we focuse on the system without spacers. 



E. Suppression of chaotic scattering contribution 

Controlling chaotic scattering trajectories requires the 
tuning of the surface disorder and the length of the mir- 
ror system. In particular, the chaotic part of the phase 
space F c can be suppressed by increasing the number of 
wall collisions as it was discussed in section HID. En- 
hancing the number of bounces is accomplished by ei- 
ther increasing the length l x or decreasing the width W. 
Furthermore, a variation of the amplitude A and corre- 
lation length l cor allows to modify the number of wall 
collisions and, moreover, to control scattering into states 
with v± > v cr . 

The sensitivity to the variation of the available con- 
trol parameters is shown in Fig. 9. While the fraction 
of chaotic phase space F c strongly decreases with l X: the 
variation with W is much weaker. Accordingly, F c does 
not simply scale with W/l x . This is, in fact, due to the 
presence of additional relevant length scales, A and l cor . 
F c and, thus, the chaotic scattering contribution outside 
of the filter band pass can be minimized for a long mir- 
ror system (large l x ) with short-ranged surface disorder 
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FIG. 9. Dependence of the fraction of chaotic scattering F c (solid line), skipping (1 — F c ) (dashed line) trajectories and the 
total transmission probability T through the neutron filter (dotted-dashed line) as a function of (a) the filter length l x \ (b) of 
the slit width W] (c) of the amplitude A and (d) of the correlation length l cor of the surface disorder. The parameters of the 
filter (if not varying): W — 250 firn. l x — l y — 10 era, A — 3 firn. l cor — 2 \ivn. 



(small I cor ) and large disorder amplitude A. Clearly, the 
price to pay for this purification of the beam is the re- 
duction of the overall transmission probability (typically 
< 10 %). 

Since a variation of the amplitude or the correlation 
length of the surface disorder might be experimentally 
difficult to achieve, a practical approach is to vary the 
length of the mirrors. In Fig. 10 we present distributions 
of the total velocity and its Cartesian components v x , v y 
and v z in the absence of the spacers and for the doubled 
length (l x = 20 cm) of the ARMS compared to the pre- 
vious settings (Fig. 6). The initial angular distributions 
entering the CTMC simulation are, as before, cosine dis- 
tributions (results for uniform distributions are hardly 
different). 

In the data shown in Fig. 10, contribution from chaotic 
scattering trajectories to the velocity distributions is sup- 
pressed and the spectra consist predominantly of skip- 
ping trajectories (see Fig. 6 for comparison). Although 
the high velocity part of the spectra (see Fig. 10 (a, b)) 
outside the upper filter limit has disappeared, there is 
still a small contribution of chaotic scattering trajecto- 
ries at velocity below the lower cut-off of the filter. This 
part corresponds to neutrons with an energy smaller than 
the Fermi potential of the mirrors and with a number 



of bounces not high enough to be absorbed due to re- 
peated wall collisions. Because of an increasing number 
of bounces the distribution within the regular island at 
the detector position, shown in Fig. 7(b), becomes more 
homogeneous compared to Fig. 7(a). From the fraction of 
the total number of detected neutrons (Fig. 9(a)), it can 
be seen that the accumulation time for good statistics is 
of the same order as for the setup previously discussed in 
section III A. 



IV. CONCLUSION 

In the present work, we propose and analyze an 
absorbing-reflecting mirror system which can be used 
as velocity band pass filter for ultracold neutrons. Full 
three-dimensional Monte Carlo simulations of the propa- 
gation of UCN through the system have been performed. 
The efficiency of the filter is explored as a function of the 
geometric control parameters and the underlying scatter- 
ing phase space structure is analyzed. 

Due to the mixed phase space, there are two types of 
motion in the system: the regular skipping motion which 
corresponds to the trajectories bouncing off the lower 
wall only and chaotic scattering trajectories which scat- 
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FIG. 10. Classical trajectory Monte Carlo simulations for the filter as in Fig. (6), however with doubled length l x = 20 cm 
and without spacers. Distributions of: (a) the absolute value of the total velocity f(v); (b-d) the velocity components f(v x ), 
f(v y ) and f(v z ) (solid line). Contributions from gravitational skipping and chaotic scattering trajectories are indicated by 
dotted blue and dashed green lines, respectively. Vertical lines correspond to analytical results for the velocity band pass (see 
Eq. (14-16)). The parameters of the filter: W = 250 /im, l y = 10 era, A = 3 fim, l cor = 2 fim. 



ter also off the upper plate with surface disorder. These 
different classes of motion give different contributions to 
the total velocity distribution of the detected neutrons. 
Since skipping motion is regular and well controlled, the 
transmission is allowed only for a limited window of the 
full vectorial velocities giving rise to a band pass filter. 
Its limits can be determined analytically. They depend 
only on the geometric parameters of the system and are 
independent of the initial velocity distributions at the 
exit of the neutron guide which are usually only partly 
known. 

The second part of the total velocity spectrum consists 
of neutrons coming from the irregular part of the phase 
space. The fraction of such neutrons strongly depends 



on the number of scattering events inside the mirror slit. 
This contribution extends over the whole range of ve- 
locities available at the collimator including those with 
high, unwanted velocities. That is why in experiments 
in which a limited velocity range is important, it is nec- 
essary to suppress the fraction of such chaotic scattering 
trajectories. We have shown that by optimizing the geo- 
metric parameters of the filter or the surface disorder the 
fraction of neutrons arriving at the detector originating 
from the chaotic sea of the phase space can be reduced 
to 5% (or less) of the total number of the detected neu- 
trons. The proposed system delivers the neutrons coming 
mostly from the regular island with a well-defined ve- 
locity range while suppressing the trajectories from the 
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irregular part of the phase space. We expect the appli- 
cation of such a filter system to UCN experiments on 
quantum scattering phenomena [22] in the near future. 
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